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ABSTRACT 

Fundamental parameters and time-evolution of mass loss are investigated for post-main-scqucnce 
stars in the Galactic globular cluster 47 Tucanae (NGC 104). This is accomplished by fitting spectral 
energy distributions (SEDs) to existing optical and infrared photometry and spectroscopy, to produce 
a true Hertzsprung-Russell diagram. We confirm the cluster's distance as d = 4611^200 PC and age as 
12 ± 1 Gyr. Horizontal branch models appear to confirm that no more RGB mass loss occurs in 47 
Tuc than in the more-metal-poor uj Centauri, though difficulties arise due to inconsistencies between 
the models. Using our SEDs, we identify those stars which exhibit infrared excess, finding excess 
only among the brightest giants: dusty mass loss begins at a luminosity of ^ 1000 L©, becoming 
ubiquitous above L — 2000 Lq. Recent claims of dust production around lower- luminosity giants 
cannot be reproduced, despite using the same archival Spitzer imagery. 

Subject headings: stars: mass-loss — circumstellar matter — infrared: stars — stars: winds, outflows 
— globular clusters: individual (NGC 104) — stars: AGB and post-AGB 



1. INTRODUCTION 

Circumstellar dust production, and the variation of its 
occurrence with a star's fundamental parameters, must 
be understood if we are to gain insight into both Galac- 
tic ecology and the history of the chemical enrichment 
of the Universe. The dusty winds of asymptotic giant 
branch (AGB) stars dominate the production of interstel- 
lar dust at all redshifts at which AGB stars are observed, 
and have therefore determined the chemical enrichment 
of Population I stars, including the Sun and Solar S ystem 
jGehrz 1989; Scdlmayr 1994; Zinncr 2003; Valiante'SaD 
[200 90 . The integrated mass loss of RGB stars determines 
the envelope mass of a star leaving the RGB tip. It conse- 
quently determines its position on the horizontal branch 
(HB) and appears to be the major factor causing the 
'second parameter' (a fter metallicity) required t o define 
HB morphology (e.g. iRool [1971 ICatelan|[2000l ). Mass 
loss on the asymptotic giant branch eventually ejects the 
star's entire hydrogen envelope, creating a post-AGB star 
and (perhaps) a planetary nebula (PN). 

A number of factors hamper the quantitative exam- 
ination of dust production. These including difficul- 
ties in determining outfiow velocities in faint, metal- 
poor stars; the difficulty in determining grain density 
(or porosity), grain size and grain shape; the diffi- 
culty in identifying the mineralogy of grains being pro- 
duced; even the difficulty in determining whether in- 
frared excess is present at all. Open questions in- 
clude how instantaneous and integrated RGB mass loss 
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varies as a function of initial mass and metallicity 
(Fusi Pecci fc Bellazzini"1997t lLamers fc Cassinellill 19991 : 
Catclan 2000) ; whether dust can form efficiently around 
RGB stars; at what evolutionary stage dust production 
begins; and whether radiation pressure on oxygen-rich 
dust grains is su fficient to dr ive a wind, especially at low 
metaUicity (e.g. lLewi^ll989| ). or whether chromospheric 
(magneto-acoustic) driving or pulsation may also be 
impo rtant (Hartmann & MacGrcgor '1980':'Du pree et al.l 
[1981 [Woitkc 2006 ; McDonald fc van Loon,2007r 

Globular clusters provide an excellent laboratory in 
which we can examine these questions. They typically 
host a single or dominant population of stars at identi- 
cal ages and metallicities, with sufficient number to be 
statistically useful. Comparisons within globular clus- 
ters therefore probe variations with evolution, while com- 
parisons between globular clusters probe variations with 
metallicity and age. To provide a proper comparison, 
it is vital to compare fundamental parameters of stars. 
These can be determined by comparing spectral energy 
distributions (SEDs) to synthetic spectra from stellar at- 
mosphere models. This approach has the advantage of 
allowing one to identify infrared excess, characteristic of 
dust production, in the SED. 

In this study, we examine the fundamental pa- 
rameters of stars in the globular cluster 47 Tu- 
canae (NGC 104), one of the most- massive (6- 
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V) = 0.04 mag (Harris 1996). T hroughout this paper, w e 
will refer to two similar studies: [McDonald et al.l (|2009D . 
which covers the cluster uj Centauri (NGC 5139; here- 
after Paper I); and [Bover et al.l (|2009| ). which covers the 
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Fig. 1. — Spatial coverage of the surveys used in this paper. See 
text for details of each survey. The cross marks the cluster centre. 

cluster NGC 362 (hereafter Paper II). 

Dust production i n 47 Tuc is of part icular interest 
due to the study of lOrigha et all ()2007[ ). This study 
has claimed observational evidence of dust production 
along the entire RGB, in contrast to the observations of 
other clusters (e.g. ui Cen; Paper I). If this claim can be 
confirmed, it significantly increases the number of stars 
we know of that are returning dust to the interstellar 
medium. This would be particularly important at high 
(z ~ 6) redshift, where AGB stars remain the dominant 
dust producers, but cannot account for a galaxy's en- 
tire dust budget (Valiantc et al. 2009). Such claims must 
therefore be taken seriously and examined carefully. An 
analysis of t he Or iglia et al. claim was undertaken by 
iBover et all (|2010[ ). who raise concern that apparently- 
red stellar colours may be a result of stellar blending in 
th e densely populated cluste r core (^-^50000 Mq pc~^ 
— iGebhardt Sz Fischerl 'lOgSl and data a rtifacts due to 
saturation of the brightest stars. Ori glia et al.l ()2010l ) 
defended their original hypothesis, suggesting that this 
dust is too warm to present considerable reddening of 
the [3. 6] -[8] colour, but does present considerable excess 
in Ks-[8]. One of our aims herein must therefore be to 
examine this issue in more depth. 

This paper follows a similar approach to Paper I. Sec- 
tion [5] details the photometric data we use, its sources 
and its reduction; determination of stellar parameters 
for the cluster's stars by fitting their spectral energy dis- 
tributions (SEDs); and the creation of a Hertzsprung- 
Russell diagram (HRD). Section [3] covers the fitting of 
stellar isochrones to the HRD. In Section 21 we derive 
mid-infrared (mid-IR) excesses for our stars, based on 
our SEDs, and use these to deter mine which stars ar e 
dusty and compare them to those in lOriglia et al.l ()2007[ ). 
The mass-loss rates and dust compositions of these stars 
are analysed in an accompanying work (Paper IV). 

2. DATA REDUCTION 
2.1. The input data 



As 47 Tu c covers a lar ge area on the sky (its tidal ra- 
dius is 43'; lHarrisl[l996[ ). a uniform set of photometric 
data is not available for every star in the cluster. The 



photometric data used in the SEDs come from a vari- 
ety of different sources covering different fields of view. 
These data are summarised in Figure [1] and Table [TJ 
They in clude data taken w ith the Spitzer Space Tele- 
scope 's (Werner et al. 2004") four Infrared Array Camera 
(IRAC; Fazio et al. 2004) bands, at 3.6, 4.5, 5.8 and 8 
fim from the Surveying the Agents of Galaxy Evolution 
in the Small Magell anic Cloud (SAGE-SMC) program 
(jCordon et al.l I2009D . At the distance of 47 Tuc, the 
1.66-1.98" full- width at half-maximum (FWHM) of the 
point spread function (PSF) of IRAC provides a resolu- 
tion of 0.036-0.043 pc, with a dynamic range sufficient 
to detect the cluster's HB stars. F urther IRAC d ata o f 
the cluster core were presented in lOriglia et al.l (|2007[ ) 
{Spitzer PID20298: PI: R. Rood]. He re, we use the re- 
analysed data from IBover et al.l ( 2010I ). For this second 
IRAC dataset, we used the "deep" photometry when [3.6] 
> 11 mag, and the "shallow" photometry otherwise. 

Longer-wavelength 24-/im data was sourced from the 
Multiband Imag ing Photo meter for Sp i tzer (MIPS; 
iRieke et al.ll2004[ i data from lBarmbv et al'l (|2009[) . as the 
SAGE-SMC MIPS data do not cover the cluster. 

Near-IR JHKs-hand photometry for the entire clus- 
ter was taken from t he 2-/im All Sky Survey (2MASS; 
iSkrutskie et al.ll2006D . The 2MASS image near the clus- 
ter core is heavily blended, and the recorded photo- 
metric fluxes are artificially raised by blending of both 
partially-resolved and unresolved stars. To solve this is- 
sue, we preferentially us ed the JHKs-hsLiid photometry 
from lSalaris et al.l ()2007| ) over 2MASS. These data cover 
a 4.'9 X 4f9 region around the cluster core. These data 
are of substantially better resolution (< 0'.'9) and con- 
sequently suffer less from source confusion and blending. 
The "deep" Salaris photometry was used in preference to 
the "shallow" when Kg > 12.5 mag. Some stars, in the 
southern half of the cluster core, were not covered by the 
Salaris data: many of these therefore have abnormally 
high 2MASS fluxes, and therefore scatter to higher lu- 
minosities. Here they can be confused with AGB stars, 
though they typically lie slightly above the AGB branch 
(their locus lies on the dividing line between regions (1) 
and (2) in the HRD shown later). 

The mid-IR Spitzer photometry is complemented by 
optical data fro m the Magellanic Cl ouds Photometric 
Survey (MCPS; iZaritskv et"aLl I2002D . which contains 
Johnson UBV-hand and Gunn i-band photometry of the 
region. Once again, this only covers the south-cast of the 
cluster. To cover the remainder of the cluster, we used 
the Jo hnson UBV -h a, iad an d Cousins /c-band photom- 
etry of iSalaris et all (|2007D, and th e Johnson UBVIj- 
bancfl photometry of lStetsonI (|2000l ). 

2.2. Creation of a master catalogue 

These data wer e com b ined using DAO- 
match/daomaster (iStetsonI 119931 ) and objects 
detected in only one fllter (assumed to be bad data) 
were dropped. The final source list contains photometry 
for 104153 stars located wi thin 30' of t he cluster core 
(00'*24"05.2" -72°04'51" — lHarrisl[T996l ). 

From this source list, a subset of sources were selected 
which had enough photometry to determine their tem- 

^ We refer to the Johnson / band as Ij throughout to avoid 
confusion with the other /-band data in this work. 
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TABLE 1 

Input photometric data and their coverage of the cluster. 



Wavelength 



Source 



Coverage 



UBVi MCPS: Zaritskv et al. (2002) 

UBVIc Stetson (2000) 

UBVIjJHKs Salaris et al. (2007) 

JHKs 2MASS; Skrutskic e t al. (2006'! 

3.6-8 /imi SAGE-SMC; (Gordon et al., in prep.) 

3.6-8 /im Origlia e t al. (2007)^ 

24 Barmbv et al. (2009) 



West of cluster, avoids core 
~ 20' X 20' around core 
Northern two-thirds of core 
Entire cluster 

South-west of cluster, not including core 
Cluster core & immediate west & east 
Cluster core & immediate south-west 



^Spitzer IRAC 
■^We use the re- 
^Spitzer MIPS; 



at 3.6, 4.5, 5.8 and 8 um: ll^'azio etaU ll200l) . 
analysed data of lBover et al.l II2010I ). 
IRieke et aO l|200l '). 



TABLE 2 

Adopted central wavelengths and zero points for the 
photometric filters used. 



Filter 


Central wavelength 


Zero point 




(nm) 


(Jy) 


U 


350 


1659 


B 


442 


4130 


V 


550 


3810 


i 


786 


2427 


Ic 


810 


2520 


Ij 


880 


2635 


J 


1240 


1602 


H 


1650 


1010 


Ks 


2160 


630 


[3.6] 


3600 


280.9 


[4.5] 


4500 


179.7 


[5.8] 


5800 


115.0 


[8] 


8000 


64.13 


[24] 


24000 


7.14 



peratures and luminosities by fitting their SEDs. We 
selected stars with at least four photometric bands with 
flux measurements, the bluest of which must be in the 
optical (/-band or shorter wavelength) and the reddest of 
which must be in the IR (iCs-band or longer wavelength) . 
This reduced the source list to 47 727 stars with usable 
photometry. 

2.3. The data reduction process 



SEDs were created following the method used in Pa- 
pers I and II. Briefly, this involves comparing broad- 
band photometry to artificially-dereddened model spec- 
tra which have been convolved to the photometric filters 
of the observations. H ere, we again use the MARCS mod- 
els of Paper I (see also lGustafsson et al.lll975ll2008l ). In- 
terpolating between these models determines the stellar 
temperature, and the multiplication of the model flux re- 
quired to match the observed flux determines the stellar 
luminosity. We initially assume that the distance and 
reddening to 47 Tuc are 4500 pc and E{B ~ V) = 0.04 
mag, respectively. We later show by isochrone fitting 
that these values are indeed appropriate (fJSl- 

We have revised the process used in Papers I and II 
to better account for int erstellar reddening. Th is follows 
the same procedure as in lMcDonald et all (120101 ) naru ely 
that we use the absorption profiles of iMcClurd ()2009l) to 
obtain A\/Aks, assuming that Ay = 3.2E{B — V) and 
that Av/Aks = 7.75. While the change in process does 
not significantly affect the stellar parameters we derive, 



it does better account for the 10-fim interstellar absorp- 
tion peak. This can affect the mass-loss rates we derive, 
though is unlikely to be significant for the small extinc- 
tion toward the cluster. 

Obtaining the correct filter transmissions is important 
when obtaining accurate stellar parameters. Differences 
can occur in the transmission efficiency of the filter it- 
self, the CCD, the telescope and (particularly) the atmo- 
sphere. We can test whether our filter transmissions and 
stellar atmosphere models are correct by plotting the ra- 
tio of the observed flux to the flux of our stellar model 
for each star (Figure [2]), which is very sensitive to errors. 
Incorrect zero points for a particular magnitude system 
manifest themselves as global offsets from unity. Incor- 
rect filter transmissions lead to a very strong tempera- 
ture dependance in the above ratio. We list the adopted 
zero points in Table [51 which gives the approximate cor- 
responding wavelengths for each filter for reference. 

Errors in one band will propagate themselves in reverse 
onto neighbouring bands, making identifying sources of 
error difficult. Such issues exist for the very hottest 
and coolest (most-luminous) stars (>6500 K and <3900 
K, respectively). At low temperatures, dynamic atmo- 
spheric processes and incomplete molecular opacities in 
our model atmospheres introduce variations from the 
model. Conversely, at high temperatures, our models are 
limited to 6500 K and do not cover the very hottest stars. 
Stars offset in Figure [2] in all bands may also be Galactic 
foreground or background SMC field stars. These are not 
well modelled as their metallicities and surface gravities 
are substantially different from those of the cluster. 

Several other systematics are visible in Figure[2] These 
include: 

• Quantisation errors among the lowest-luminosity stars, 
which manifest themselves as fan-shaped spreads visi- 
ble in the U and 8-/zm filters, but are negligible among 
the hottest stars. 

• Poor convergence of models around 6250 K. This only 
affects the hotter stars which have no mid-IR photom- 
etry, for which we cannot therefore determine mid-IR 
excess. 

• Increased scatter in the HB stars around 50 Lq. This 
is due to unknown processes (though could indicate 
an abundance spread) and does not affect any giant 
branch dust producers. 

• Small (~1%) remaining systematic under- or over- 
estimates of fiux in several filters. The colour-based 
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Fig. 2. — Differences between observed and SED-fitted model fluxes as a function of luminosity (top panels) and temperature (bottom 
panels, where only stars with L > 5 Lq are shown). 
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TABLE 3 

Stellar parameters as determined by SED fitting. A full 

ELECTRONIC TABLE IS AVAILABLE ONLINE. 



Sequence 


RA (dog) 


Dec (deg) 


Temperature 


Luminosity 


Number 


{J2000) 


(J2000) 


(K) 


(Le) 


1 


4.641571 


-72.328035 


6147 


258.4 


2 


4.685675 


-72.365643 


4449 


56.53 


3 


4.706252 


-72.387843 


5399 


1.625 


4 


4.722603 


-72.374671 


4500 


1.228 


5 


4.725042 


-72.366471 


5308 


1.057 



detection criterion we use are chosen to maximise 
detection of IR excess in such cases, but the amph- 
tude of these effects is much smaller than the >0.1 
mag (>10%) excesses we are looking for. 

We therefore do not expect these systematics to affect 
our detection of circumstellar dust. 

Problems caused by unresolved blending in the cluster 
core are visible in several bands. As noted in the previous 
section, this is a particular concern in the 2MASS data, 
where there is considerably more scatter in the fitted 
parameters and observed: expected flux ratios. 

Table [3] contains the list of stellar parameters deter- 
mined for the stars. 



2.4. Accuracy of the fits 
2.4.1. Analysis of departures from the models 



As noted in ^2.'S[ the goodness-of-fit of our SED mod- 
els depends strongly on the zero points and filter trans- 
missions used. While Figure [2] shows that the models 
agree well with the observations between temperatures 
of 3900 and 6500 K, stars lying outside of this range 
are subject to some error. To assess the magnitude of 
this error, a representative selection of stars was taken 
and the data from the most-affected filters {U, B and 
Ij) were removed. This had a negligible (<5 K) effect 
on the stellar temperatures of the stars in the temper- 
ature range 3900-6500 K, with temperatures changing 
slightly more for stars outside this range. The coolest 
stars have temperature uncertainties caused by poorly- 
matching models of order tens of degrees, but these stars 
are also strongly variable. Any 'instantaneous' measure 
of temperature such as this will therefore not necessarily 
be representative of the mean stellar temperature to this 
degree of accuracy. 

A formal estimate of the error may be calculated using 
'simulated' observations using the Monte-Carlo method, 
whereby the fiux in each band is set to be the observed 
value plus-or-minus an error sampled from a Gaussian 
probability distribution with a identical to the reported 
photometric error. This suffers two problems: (1) it 
would take prohibitively long to perform this test for 
every star and (2) it only takes into account the photo- 
metric error, not the modelling error. We have therefore 
selected a small sample of stars from across the HRD, 
and performed Monte-Carlo simulations to determine the 
errors in their parameters. We do this twice: the first ac- 
counts for only the photometric error, as described above; 
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Fig. 3. — Monte-Carlo fitting results, showing (left) photomet- 
ric errors only and (right) errors derived from departures from the 
model. Red bars show the 1-D standard deviations of the distri- 
butions. Labels denote sequence numbers of the objects (Table [Sjl, 
chosen to be representative of their respective parts of the HRD. 
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TABLE 4 

Stars chosen for Monte-Carlo simulations. 









' 1 V>ryi noval"! IT'O"'" 


T ,1 1 m 1 nrit;! 

Li-lllIllvJol L y 


Coverage 




Number 


(J2000) 


(J2000) 


(K) 


(Lq) 








12515 


5.9604583 


-72.0725000 


3916±4 ±25 


982±11 ±35 


B - 


[24] 


Upper RGB/AGB 


9352 


5.9282500 


-72.1920278 


4265±5 ±55 


316±3 ±7.2 


B - 


Ks 


Central RGB/AGB 


4929 


5.7577500 


-72.1565278 


4567±6 ±59 


184±1.4 ±16 


B - 


[8] 


Lower AGB 


8937 


5.9237083 


-72.0840556 


4329±10±65 


175±3 ±11 


B - 


[24] 


RGB 


19514 


6.0277917 


-72.0754167 


5283±30±80 


60.0±0.7 ±1.7 


B - 


Ks 


HB (core) 


28472 


6.1138333 


-72.0774722 


5191±27±18 


59.8±1.0 ±0.5 


B - 


Ks 


HB (not core) 


15423 


5.9881250 


-71.9649722 


4782±10±36 


45.0±0.4 ±0.6 


B - 


Ks 


RGB clump 


40122 


6.5137500 


-72.2503611 


6250±13±56 


5.05±0.06±0.18 


B - 


[4.5] 


MSTO 



^ First error encompasses photometric errors, second error encompasses modelling errors. 



the second replaces the Gaussian a with the deviation of 
the observations from the model. 

The stars with which we perform this test are listed in 
Tableland the results are shown in Figure [31 It is firstly 
clear that the temperature and luminosity we derive are 
interdependent: spread along the RGB track appears de- 
termined primarily by errors in the optical data, while 
spread perpendicular to the RGB track appears deter- 
mined primarily by errors in the IR data. It is also clear 
that the photometric errors are considerably less than the 
total error as estimated by differences between the ob- 
servations and stellar atmosphere models. As we note in 
^4.11 this is probably due to under-reported errors due 
to image artifacts and blending from sources within a 
few arcseconds. In the right-hand panels of Figure [31 we 
note that those stars with no data longward of if^-band 
have significantly more scatter perpendicular to the giant 
branch: Spitzer data are therefore crucial in determining 
accurate stellar parameters, even for warm stars. 

Varying a star's fundamental parameters can also have 
an effect on the values we achieve. Likely variations for 
bona fide 47 Tuc cluster members are somewhat smaller 
than the corresponding errors in temperature or lumi- 
nosity caused by inaccuracies in the model, which are 
of order 1-2% for temperature and 2-9% in luminosity. 
By varying [Fe/H] between solar and -1.4 we produced 
differences of similar magnitude, varying the mass from 
0.6-1.0 Mq produced differences of order 0.3% in lumi- 
nosity, and varying E{B — V) between and 0.12 mag 
produced luminosity changes of <5%. The largest un- 
certainty, at least in terms of luminosity, appears to be 
in the distance to the cluster. 

2.4.2. Comparison with other works 

In order to determine the accuracy of our SED 
fits, we compare our temperatures with a number of 
spectroscopically-derived temperatures. We do not com- 
pare the luminosities, as the fitting procedure means that 
the luminosity is directly determined from the model af- 
ter the temperature fit has been made. The resulting 
discrepancies are listed below: 



• iCarretta etHI (|2004[ ) — nine early-RGB stars («5100 
K), temperatures agree within errors: on average they 
are 21 K (0.4%) warmer in our SED fits than their 
spectroscopic temperatures, cf. 108 K for the standard 
deviation of the temperature differences. 



atures also agree within errors, averaging 47 K (1.3%) 
cooler in our SED fits (st. dev. 311 K). 



i [McDonald fc van LoonI (|2007l ) — nme upper 
RGB/AGB stars (<4000 K, >1100 Lq), temper- 



• iKoch fc McWilliamI (|2008| ) — eight central RGB stars 
(4200-4500 K), individual temperatures differ at the 
2(7 level. As a whole, the stars average 25 K (0.6%) 
warmer in our SED fits, smaller than the standard 
deviation of the differences between the two papers 
(31 K). 

• 'Worlcv et al.l (120101) — five of the spectra used in 
iMcDonald fc van LoonI (|2007[ ). individual tempera- 
tures differ at the 4(t level, averaging 231 K (6.1%) 
cooler in our SED fits (st. dev. 108 K). 

The exceptionall y low standard deviation of the 
iKoch fc McWilliani (|2008D temperatures provides a good 
check of the accuracy of both our measurements and 
theirs. It is worth noti ng that the discrepan cy between 
our results and those of ICarretta et al.l (|2004l ) lies in the 
accuracy of our photometry, not their spectroscopic de- 
terminations. A systematic offset between our SED fit- 
ting and conventional spectroscopic techniques is likely 
hmited to <30 K (<0.7%) in temperature and <2.8% in 
luminosity. 

2.5. The Hertzsprung-Russell diagram 

Having determined the stellar parameters for the above 
objects, the data were then compiled into an HRD, pre- 
sented in Figure m 

The HRD shows a number of features, mostly iden- 
tifying our target stars as post-main-sequence stars in 
47 Tuc. Also visible on the HRD is a background se- 
quence of stars belonging to the SMC (13). This merges 
with the main sequence (MS) of 47 Tuc at roughly the 
luminosity where the SMC's HB and red clump would 
be. Regions (9) and (14) contain foreground and back- 
ground objects of indeterminate distance. Typically, we 
would expect region (9) to contain stars in the Galactic 
disc (47 Tuc li es at Galacti c co-ordinates I — 305.90°, 
b = -44.89° — lHarrislll996D and regio n (14) to contain 
background galaxies ([Bover et al.ll2008t Paper I). In this 
particular sample, however, region (14) may also contain 
some SMC young stellar objects (YSOs) and cool main 
sequence objects. We do not expect many YSOs, how- 
ever, as our region of interest lies away from the bulk 
of the galaxy's star-forming regions. The main sequence 
stars at this temperature are mostly below the sensitivity 
threshold of our data. 

Vertical artifacts are also visible on the diagram. These 
occur every 250 K and correspond to the temperature 
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Fig. 4. — HRD of stars in the direction of 47 Tuc, showing evolutionary and spatial components. The top panel contains all stars, the 
bottom panel contains only stars more than 2' from the cluster core, which have photometric data in the optical, near-IR and mid-IR. 
Abbreviations are discussed in the text. 



8 



McDonald et al. 



grid of our model spectra. They are generally caused by 
poor-quality optical photometry in low-luminosity stars. 
Especially important are discrepancies between the ob- 
served U- and i?-band photometry and those predicted 
by the MARCS model spectra, hence the increased preva- 
lence of artifacts at higher temperatures where more flux 
is emitted in U and B. As our models only extend to 6500 
K, we must assume that stars with temperatures derived 
to be >65GG K have progressively more-uncertain param- 
eters. We have therefore limited our study to stars with 
T < 7000 K, of which there are 46 398. 

Several important evolutionary points (numbered in 
Figure^) can be noted at this point, which we use in the 
next section to fit stellar isochrone models. These are: 
the main-sequence turnoff (MSTO; 11 in Figure [J]), the 
start of the RGB (10), the RGB bump (8), the upper 
RGB (5), the RGB tip (top of (2)), the base of the HB 
(bottom of (7)) and the start of the AGB (bottom of (4)). 
A review of these evolutionary processes, wit h particular 
emphasis on globular clusters, can be found in lMcDonaldl 
(|2009). 

Circumstellar dust reprocesses optical light into the IR, 
causing an apparent cooling of the star's effective temper- 
ature. This moves the star to the right on the HRD, mov- 
ing into region (2). Throughout their time on the RGB, 
stars lose mass through stellar winds. For the majority of 
the time, these appear to be dus tless winds driven from 
the chromosphere (M cDonald fc van Loon 2007; Paper 
I), but if the IR excesses reported by Origlia ct al. (2G03) 
are true, we may see some low-luminosity stars in this re- 
gion by virtue of their circumstellar dust. This issue is 
confused, however, by the presence of field stars, and of 
stars with poorly-determined photometry. 

Our HRD contains a well-denoted and well-separated 
early-AGB. By its very nature, however, the AGB 
asymptotically approaches the RGB, and the two become 
indistinguishable above ~300 Lq. More-rapid evolution 
on the AGB means that there are 4-5 times more stars 
per unit luminosity on the RGB than AGB. 

Helium fusion on the upper AGB also becomes increas- 
ingly volatile: its violence can increase a star's luminosity 
dramatically (up to a factor of several) for a short period 
of t ime (10^ ^ 10^ years), termed a th ermal pulse (TP; 
e.g. IZiilstral [19951 : IWachter et al.l[2008h . Such stars may 
scatter above the main giant branch in the HRD, into re- 
gion (1), but this area is confused by normal AGB stars 
with poor photometry. 

Stars on the TP- A GB are also unstable to pulsation 
via the k- mechanism (jUlmschneided IT998I ) . where har- 
monic oscillations are stochastically excited. This causes 
semi-regular variability at first, then more-regular and 
increasingly-violent pulsation until the AGB tip, where 
the stellar atmosphere disperses. This variability will 
cause some horizontal and vertical scatter of these stars 
in the HRD, as we typically only use one epoch of pho- 
tometry for every filter. 

3. ISOCHRONE FITTING 

3.1. Introduction 

By fitting stellar isochrones to our data, we can re- 
determine the basic parameters for the cluster as a whole, 
namely its distance, reddening and age. To do this, we 
must make assumptions about the cluster's metallicity. 
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Temperature (K) 

Fig. 5. — Top panel: Padova isochrones generated using the CMD 
utility, overlaid on a density-mapped HRE). Bottom panel: Padova 
isochrones generated using the YZVAR utility. Arrows show the 
locations of the RGB bumps of the isochrones. Details of the 
isochrones' parameters can he found in ^3.11 

helium fraction and a-element enhancement. The accu- 
racy of this determination depends both on the system- 
atic errors present in our data fi 32.4|) . and the validity of 
the assumptions we make. To minimise systematic er- 
rors, we have used the 'cleaned' subset of stars presented 
in the bottom panel of Figure 21 to the exclusion of stars 
with poor photometry and those in the immediate cluster 
core. 

In this study, we compare our HRD to two sets of 
isochrones, both of which were use d in Paper I. These 
are n amely th e Padova isochrones (jMarigo 
Bert elli et al.l [2008j3 and the Dartmouth isochrones 
(Potter et al .l[2008fl. 

The CMD inte rface for the Padova isochrones 
(jMarigo et al.l l200"8l ) does not allow the user to change 
either the helium content or g -enrichment of th e 
isochrones. The YZVAR interface (jBertelh et al.l 12008V 
however, allows bot h the helium co ntent and a Reimers- 
law mass-loss rate (jReimersI I1975D to be varied by the 
user, though also does not allow the a-enrichment to be 
varied. In this regard, we have simply adopted a fixed 
metal content, [Z/H], and assumed that metals scale as 
solar values. The Dartmouth models allow the variation 
of both helium and a-enrichment in fixed steps. 

Altering the above parameters has the following effects: 

• Increasing the helium content of the isochrones steep- 

^ http:/ /stev. oapd.inaf.it/cgi-bin/cmd 

^ http:/ /stellar. dartmouth.edu/~models/webtools. html 
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Fig. 6. — As Figure[5]for the Dartmouth isochrones. Top paneh 
isochrones at [a/Fe] = 0. Bottom paneh isochrones at [a/Fe] = 
+0.2. 



ens the slope of the RGB, though the RGB tip stays in 
roughly the same place. It also causes an increase in 
the luminosity difference between the MSTO and the 
start of the RGB. 

• Increasing the a-enrichment accelerates MSTO, mean- 
ing this point corresponds to lower luminosities. It also 
shifts the entire isochrone to cooler temperatures and 
increases the luminosity of the RGB bump. 

• Increasing the metallicity of the isochrones will shift 
them to cooler temperatures, as stronger photospheric 
lines make the atmosphere more opaque. 

• Increasing the fitted distance of the cluster will shift 
the isochrones toward higher luminosities with respect 
to the stars. 

• If the cluster's reddening is higher, we will have there- 
fore underestimated the de-reddening correction: fur- 
ther correction then requires the stars to move to 
higher temperatures and luminosities with respect to 
the isochrones. This does not affect all stars equally, 
and hotter stars will require larger corrections than 
cooler stars, with the result that the HRD becomes 
horizontally stretched and vertically skewed (see Pa- 
per I for the way this is treated). 

• The cluster's age affects the isochrones in a number 
of ways, with the largest alterations being in the posi- 
tioning of the MSTO and base of the RGB. 

We initially assumed the following parameters: 



• MetaUicity: [Z/H] = [Fe/H] = -0.7; 

• Hehum fraction: Y = 24% « Yq (jPorman et al.lll989[) : 

• a-enrichment: [a/Fe] ~ +0.2 dex (Dartmouth 
isochrones only); 

• Distance: d — 4.5 kpc; 

• Reddening: E{B - V) = 0.04 mag; Ay = 0.124 mag; 
R = 3.1; 

• Age: t = 12.5 ±li Gyr. 

Isochrones using these parameters are shown in Figures 
[5] and HI It is clear from these figures that the best- 
fitting isochrone is that of the Dartmouth models, but 
at the adopted value of [a/Fe] = 0. This provides a 
remarkably-good fit to the data. The only exception to 
this is on the upper RGB, above the red clump, where 
the [a/Fe] = +0.2 model fits better (see also later, in 
Figure [9]). An increase in a-element abundances follow- 
ing first dredge-up is not expected, as the stars have not 
begun helium burning. It could, however, be explained 
if extra mixing on the RGB occurs around first dredge- 
up, and if such mixiiig alters the surface chemi stry (e.g. 
iD'Antona fc Ventural 120071 : iKarakas et aIll2010L and ref- 
erences therein). This effect is shown in more detail in 
t j3.2.5l Note that the Padova isochrones, which qualita- 
tively fit the upper giant branches and MSTO, are also 
for [a/Fe] 0. 

3.2. Cluster parameters 

3.2.1. Metallicity and abundances 

The metallicity and abundance of 47 Tuc's stars is bet- 
ter determined from spectroscopic studies of individual 
stars than from our HRD (see §2.4.2l for such studies). As 
noted, however, the best- fitting isochrone requires [a/Fe] 
= 0. To fit the Padova isochrones by changing the a- 
enhancement (which roughly corresponds to moving the 
isochrone horizontally in the HRD) requires [a/Fe] < 0. 
For the purposes of isochrone fitting, we have proceeded 
under the assumption that there is no a-element enrich- 
ment compared to solar values. 

3.2.2. Reddening 

The de-reddening of the cluster we apply can be al- 
tered using the relations listed in Paper I. This gives an 
approximate solution for different values of E{B — V) 
without the need to re-analyse the entire dataset. Based 
again on the Dartmouth [a/Fe] isochrones, we find agree- 
ment with th e accepted value of E{B — V) 0.04 mag 
(IHarrislll996l) . 

3.2.3. Age and distance 

Once the metallicity, abundances and reddening of the 
stars have been set, the age can primarily be determined 
from the temperature of the MSTO. This is vertically 
scaled to match the luminosity of the MSTO /RGB base 
at 5500 K. The parameters of age and distance are thus 
anti-correlated. 

As shown in Figure [71 isochrones for 12 ± 1 Gyr fit the 
data well. An incomplete stellar catalogue of objects be- 
low ~1 Lq means that the models will fit progressively 
less well below this luminosity. Isochrones are plotted 
for distances of 4783, 4611 and 4456 pc for 11, 12 and 13 
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Fig. 7. — As Figure [6l showing Dartmouth isochrones at [a/Fe] 
= 0, covering the MSTO. Isochrones are placed at 11, 12 and 13 
Gyr. 
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Luminosity (Lg) at 4.5 Wpc 

Fig. 8. — The luminosity function for 47 Tuc (solid, red line), 
along with theoretical Pisa isochrones. Lower, short-dashed line 
(blue): AGB and HB stars only. Upper, long-dashed line (green): 
AGB/HB -I- RGB stars. 

Gyr, respectively. Adding in the 2.8% systematic lumi- 
nosity uncertainty in luminosity (i |2.4.2p . this equates to 
a distance of 4611 tjoo P'^- 

Both sets (CMD & YZVAR) of Padova isochrones yield 
similar constraints on age. The distances implied for 
these isochrones are approximately 4% (180 pc) less than 
for the Dartmouth isochrones. 

3.2.4. Luminosity function 

Figure [8] shows the observed and theoretical luminosity 
functions of 47 Tuc. The observed luminosity function 
is derived from all stars below 6000 K and 4200 L0 (to 
include the dusty long-period variable star 47 Tuc VI, 
but exclude known field stars) which do not fall in region 
13 (SMC) of the HRD. The theoretical isochrone uses the 
Pisa stellar evolution models for a 0.90 Mq RGB and a 
0.65 Mq AGB star. The isochrones have been scaled to 
one star for every 80 000 years of evolution. It should be 
noted that the AGB models terminate at 1550 Lq: we 
have extrapolated them in Figure|S]to estimate that there 
are some 8±2 AGB stars above the theoretical RGB tip 
at 2780 Lq. 

Overall, the luminosity function provides a good fit to 
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Fig. 9.— The upper RGB, HB and AGB of 47 Tuc (small, black 
dots). These are overlain with Dartmouth RGB isochrones (blue 
dashed lines) at [o/Fe] = (left) and -fO.2 (right), and with similar 
HB models (large, coloured dots). In order of increasing luminosity, 
the HB tracks are for stellar masses of 0.85, 0.90 and 0.95 Mq. 
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Fig. 10. — As Fig. [21 but with Pisa HB models at (left-to-right) 
0.65, 0.70 and 0.75 Mq. 



the data. The departure between 60 and 100 Lq is prob- 
ably due to contamination from field stars and scatter in 
the HRD from poor photometry. The fact that we do not 
find as many HB stars at ~50 Lq as the models predict 
suggest that either we do not have a complete sample of 
all HB stars, or that some of the HB stars have already 
been ejected from the cluster due to mass segregation 
(cf. HB lifetime » 190 Myr; core relaxation time « 91 
Myr; lHarrisl[l996[ ). This is not readily visible in radial 
density plots of our sources due to the differing depths 
of the original surveys we used. 

3.2.5. HB mass via isochrone fitting 

The amount of mass lost on the RGB can be deter- 
mined by measuring the difference between the RGB-tip 
stars' initial masses and the mass of stars on the HB. The 
age of the stars means we can negle ct mass loss during 
the c omparatively short (~1.8 Mvr: [SiTva Aguirre et al.l 
l2010f ) evolutionary period between the RGB-tip and 
zero-age HB. Initial masses for isochrones at 12 ± 1 Gyr 
are: 



Fundamental parameters and dust production in 47 Tuc 



11 



150 



°- 100 

in 90 

^ 80 

31 70 

^ 60 



50 



40 



1 1 


1 1 1. . If;li%>» .. 

■jS'.«<V'. 




1 Iv 




.fi "i*.'.';"" ■ ■ 




A' s • 






- ■ ' 




- 


■ ^-^'"^''''^ /''•''' ■' W'- ' 








-;;^.::;-'-'f?^^-. O.75;:M0 •'^-i-- _ 




















■ 1 ■ J ■ 


1 1 .1.. ■ : 



10000 



9000 



Fig. 



8000 7000 
Temperature (K) 
11. — As Fig. 1101 for the cluster Cen. 



6000 



5000 



0.893 
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-0 021 ■'^0 for the Padova isochrones; 
at [a/Fe] = 0; 

• Minit = 0.889 10022 Mq for the Dartmouth isochrones 
at [a/Fe] = +0.2; 

• -^finit = 0.916 M(7) for th e Pisa evolutionary models 
(from iGratton et all [20ll . 

HB star masses can be estimated from synthetic HB 
models. The temperature spread of the HB is greater 
than that of the RGB at the same luminosity (the stan- 
dard deviations from the mean temperature are « 2.2% 
vs. ~ 1.6%, respectively). This implies a moderate varia- 
tion of stellar parameters among the HB stars. Notably, 
the HB lies at roughly constant luminosity, though lu- 
minosity may rise very slightly towards higher tempera- 
tures. 

Figures [51fc [TU|show Dartmo uth and Pisa synthetic HB 
models (Dott eTet al.l 1200^^1 : |Castcllani ct al. 200^3) 
for different stellar masses. The Dartmouth models 
are given at [a/Fe] = and -1-0.2; the Pisa models at 
Y — 0.24. Two things are immediately obvious from 
these figures. Firstly, that the temperature spread of 
the HB is not due to different stellar masses: the near- 
constant-luminosity morphology is more consistent with 
a spread in a-element enhancement. Secondly, the HB 
star masses predicted by the two synthetic models varies 
considerably. Including the 2.6% systematic luminosity 
uncertainty and 4.5% distance uncertainty, we estimate 
from Figures [9l fc [TOl that the mass of the HB stars is 0.90 
± 0.05 and 0.92 ± 0.05 Mq for the Dartmouth models, 
assuming [a/Fe] « and +0.2, respectively, and 0.68 ± 
0.03 M0 using the Pisa HB tracks. 

This shows that the two HB models are clearly mutu- 
ally incompatible, and that significant work in this area 
must be done before an agreement can be reached. The 
Dartmouth model suggests that RGB mass loss must be 
remarkably low, at < 0.06 (at the +l-cr level, i.e. 
for 84% of the stars). The Pisa models, however, predict 
RGB mass loss to be 0.24 ± 0.04 M© at the ±l-cr level 

http://stellar.dartmouth.edu/ models/ 
" http://astro.df.unipi.it/SAA/PEL/ZO.html 



in ag reement with the 0.242 Mq found bv lGratton et al.l 

2oToh ). 

The small range of modelled HB-star masses provided 
by each set of models also suggests that mass loss on the 
HB is low (< 0.1 Mq) in 47 Tuc, which would indicate 
that chromospherically-driven win ds are not so impor- 
tant in this situation as implied bv lDupree et al.l ((2009): 
more consistency between the HB models would be nec- 
essary to substantiate this statement. 

The integrated mass loss and median mass predicted 
by the Pisa models (saO.68 Mq) is in good agreement 
with that of the stars in w Cen (Figure [Til data from 
Paper I). The majority of w Gen's stars are 8-9 times 
more metal-poor than those of 47 Tuc, but crucially are 
of very similar age (hence similar initial masses). Due 
to the lack of consensus between the models, we use the 
Pisa models as a rough upper limit to the integrated RGB 
mass loss occurring in 47 Tuc. 

It should also be possible to determine the mass of 
HB stars directly, by taking our photometrically-derived 
luminosities and temperatures, and comparing them to 
literature spectroscopically-derived surface gravities (g), 
via: 

L GM 
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where the symbols take on their usual meanings. Due to 
the systematic uncertainties (in particular in measuring 
surface gravity), this can only be done with confidence 
on a set of stars reduced using identical processes, and 
preferably with similar temperature and luminosities. 
Unfortunately, there exist no spectroscopically-derived 
measurements of surface gravity measurements for a se- 
lection of both RGB and AGB stars with the internal 
accuracy required to do this. 

If one assumes an integrated RGB mass-loss rate of 
0.22 Mq per star, comparable with the 0.20-0.25 Mq 
derived for uj Centauri in Paper I, an AGB star will typ- 
ically have a mass of ~0.67 Mq. This implies that only 
~0.15 Mq of its envelope mass remains, to be subse- 
quently lost towards (or at) the AGB tip. We stress, 
however, that systematic uncertainties in the HB models 
mean that this cannot be proven from our data. We re- 
turn to this subject of integrated RGB mass loss in the 
discussion. 

4. MASS LOSS AND DUST PRODUCTION 

Mass loss during the TP-AGB can give rise to dust 
production in the stellar outflow. The chemistry of the 
dust is determined by the amount of dredge-up which 
occurs: this defines whether carbon or oxygen is more 
numerous in the stellar atmosphere. Carbon and oxygen 
combine in the atmosphere to form CO. The remain- 
ing carbon or oxygen then goes on to define whether 
the dust is carbon-rich (primarily in the form of amor- 
phous carbon) or oxygen-rich (primarily in the form of 
silicates). While a few carbon stars a re found in globu- 
lar clusters (e.g. Ivan Loon et al.ll2007f ). they tend to be 
rare: their low mass means third dredge-up is not usu- 
ally sufficient to cause the C/0 ratio to exceed unity, 
thus most stars produce silicate dust. Metallic iro n dust 
can also form around oxygen-rich stars ( Kemper et al.l 
[200llVerhoelst et al.|[2009t lMcDonaideral.ii201Q) . It is 
not clear whether metallic iron dust forms in carbon-rich 
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stars too, due to the difficulty of spectroscopically sepa- 
rating metallic iron from amorphous carbon dust. 

The dust can then drive a dusty wind by absorbing 
momentum from stellar radiation incident on it. This 
radiation pressure forces the dust from the star. Dust 
grains are collisionally-coupled to some extent with the 
surrounding gas, meaning that the gas is driven away 
from the star as well. It may be that stellar pulsations 
provide a net outward velocity to the dust, either as an 
initial 'kick' velocity, or by the dissipation of acoustic 
energy in the extended atmosphe re (e.g. |Bpwcn 1981; 
iLewisllTQSl Ivan Loon et"alll2008b[ ). 

Dust will re-radiate absorbed optical stellar light in 
the IR. This has the apparent effect of cooling the star 
and giving it excess in the IR above our model spectrum. 
Dust emits a modified blackbody spectrum which exists 
in addition to the stellar photospheric output at those 
wavelengths. Amorphous carbon, graphite and iron dust 
have no IR spectral features: such dust will typically give 
rise to positive values for the colours [5.8]-[8], [4.5]-[5.8] 
and even [3.6]-[4.5], though the latter depends on the 
dust's temperature. Particularly warm dust may also 
produce {Ks-[i-6]) > in addition to the above colours. 
Silicate dust leads to broad emission at 9.5 and 18 /im, 
which is partly covered by Spitzer's 8- and 24- /im filters, 
giving rise to an excess in these bands. 

4.1. The infrared excess stars 
4.1.1. Determining infrared excess 

Stars will exhibit IR excess for a variety of reasons not 
attributable to circumstellar dust. These include: 

• Increased photometric and astrometric uncertainty 
due to problems with source separation in the dense 
cluster core; 

• Photometric errors due to 'ghost' stars and linear arti- 
facts near bright stars (caused by the bandwidth effect 
and banding, respectiveljQ); 

• Artificial brightening due to blending with unresolved 
objects in mid-IR observations, which are at lower res- 
olution than their near-IR/optical counterpart obser- 
vations; and 

• Artificial brightening due to blending with other IR- 
bright objects, such as background galaxies or dusty 
SMC stars. 

Of these, the first two will increase roughly as the square 
of the stellar density, while the latter two will increase 
proportionally with the stellar density. 

To ensure we only select those stars that harbour cir- 
cumstellar dust, we wish to perform a visual inspection 
of the photometry and imagery of each star. Such an 
inspection is performed to mitigate against data arti- 
facts that may otherwise go undetected. It also allows 
us to sensibly include photometry from other literature 
sources (e.g. the A KARI IRC point source catalogue; 
llshihara et al.|[20Tol) which improve coverage of dust fea- 
tures in the SED, but where coarser resolution and/or 
lower signal-to-noise mean flux measurements need to be 
interpreted with caution on an individual basis. 

Details can be foun d in the I RAC handbook: 
|http: //ssc. spitzer.caltech.edu/irac/iracinstrumenth andbook/ 



It would take prohibitively long to do this for our orig- 
inal list of 46 398 stars. We reduce our source list to a 
manageable number by selecting only stars with L > 31.6 
Lq (i.e. the 3299 stars in or above the red clump) which 
have Spitzer photometry (1993 of those stars) and meet 
one or more of the following criteria: 

• is an evolved star candidate ( Wl.bl 58 stars); 

• has at least two of the (iC^-[3.6]), ([3.6]-[5.8]), ([4.5]- 
[8]) or ([8]-[24]) colours which are more than 0.1 mag 
greater than that predicted by the model (74 stars); 

• has {Ks-[8]) or ([3.6]-[8]) at least 0.1 mag greater than 
that predicted by the model, if more-luminous than 60 
L0, or 0.2 mag greater otherwise (129 and 80 stars for 
each colour, respectively); 

• has (iirs-[24]) or ([3.6]-[24]) excess under the above cri- 
teria (36 and 42 stars for each colour, respectively); 

• has photometry from AKARI (35 stars); 

• is a known variable (43 stars). 

When combined, this list gives 258 unique objects. One 
may wonder why we use colours relative to our model 
in this case, rather than the absolute colours. In cooler 
stars, CO bands at 2.3 and 4.5 /im cause substantial 
deficits in the Kg and 4.5-/im bands. In cool stars, 
around 3500 K, we model the colours of a naked star 
to be {K,-[^.6]) = 0.17, ([3.6]-[5.8]) = 0.05, ([4.5]-[8]) = 
0.3, {K,-[8]) = 0.33 and ([3.6]-[8]) = 0.16 mag. These 
colours decrease in warmer stars. By taking into account 
the expected colours and their variation with stellar pa- 
rameters, we can probe smaller excesses than possible 
with conventional colour-magnitude diagrams. The stars 
we investigate are shown in such a colour-magnitude di- 
agram in Figure 1121 

A large number of stars were flagged where the pho- 
tometry agreed between the 3.6- and 5.8-/im bands and 
between the 4.5- and 8-/im bands, but were discrepant 
between these two pairs. This is likely due to the si- 
multaneity of the 3.6- and 5.8-/im, and 4.5- and 8-/im 
observations. It could either be a direct consequence (in- 
trinsic variability) or indirect consequence (such as dif- 
fering qualities of photometry due to different coverage 
in the original mosaics). For this reason, we have ignored 
the ([3.6]-[4.5]), ([4.5]-[5.8]) and ([5.8]-[8]) colours in the 
above criteria. 



At this point, we have a list of potentially dusty ob- 
jects, using criteria w hich ar e qual itatively similar to 
those applied in Origli a~et al.l (|2007l ). Astrometry and 
photometry of the sources in that work have been made 
available to us (L. Origlia, private communication). Of 
their 93 sources, only 45 are found to have counter- 
parts with IR excess in our own candidate list (Table 
[S]). This immediately suggests that the primary differ- 
ence between the two works lies in the original photo- 
metric reduction, and that the real photometric errors 
ar e considerably large r than those used for the 3cr cutoff 
in lOriglia et al.l ()2007D . Of these 45 matched sources, 23 
are over 1000 L0, while only 22 are below 1000 L0. We 
explicitly do not claim one reduction to be more accu- 
rate than the other: while several stars known to harbour 
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Fig. 12. — Absolute (Ks — [8]) and {[3.6] — [8]) colours for all stars in our sample (small, red dots). Circled dots (green) show our initial 
cut for dusty candidates. Triangles (blue) show our final cut for stars which we scrutinise individually, once bad data have been cleaned. 



dust do not appear in Origlia et al.'s list (e.g. V 13 and 
V18; Ivan Loon et all 120061 : iLebzelter et all [20061) . they 
lie outside their survey region. 

One comparison we can make is of the magnitudes 
of Origlia et al.'s dusty stars used in their work, and 
those used in this work (Figure IT5)) . It is notable that 
(particularly for the fainter stars) the X^-band magni- 
tudes are systematically brighter in this work, especially 
where lower-resolution 2MASS data are available, while 
the [3.6] and [8] magnitudes are systematically fainter. 
This is further evidence that the method of photomet- 
ric reduction and d ifferences in r e soluti on a re the major 
differe nces between [Origlia et al.l (|2007[ ) and lBover et al.l 
([2OT0h . 

4.1.2. Sample cleaning 

Most of the 258 objects that match our criteria above 
were flagged as a result of incorrectly-matched photom- 
etry, either between surveys in our cross-matching rou- 
tines, or within the individual surveys themselves. The 
Gunn i-band data from the MCPS and the JHKs data 
from 2MASS are particularly prone to these errors. We 
expect dusty stars to have a mid-IR excess, and that 
the ratio of observed/modelled flux should not apprecia- 
bly decrease toward longer wavelengths. To remove stars 
which have been erroneously selected due to bad photo- 
metric data or matching, we have removed stars selected 



with the following cuts: 

• To remove objects with poor JHKs-Spitzer matching, 
we remove objects which show i^s-band excesses (com- 
pared to our model) more than 0.1 mag greater than 
the associated 3.6-/im excess. An example is shown 
in Figure [T4l (i). This reduces our source list to 221 
objects. 

• To remove objects with poor JHKs matching in gen- 
eral, we remove objects which have observed/model 
flux ratios at J, H and /C^-band of <0.8 (i.e. a strong 
near-lR flux deficit; Figure [HI (ii)). This leaves 188 
objects. 

• To remove objects with bad photometry in Spitzer 
data, while retaining stars with only 8- or 24-/im ex- 
cess, we remove objects which show less than 10% ex- 
cess above the model at the longest recorded wave- 
length: 12 sources are discounted because they show 
no excess at 24 /im and 37 because they show no ex- 
cess at 8 jjLva (when 24 and 8 ^m are the longest wave- 
length with data, respectively). These stars mostly 
have discordant IRAC photometry, but include V16, 
LWl, LW2 and LW21, which were selected due to their 
variability, but show no apparent mid-IR excess (Fig- 
ure[T4l (iii)). This leaves 128 objects. 

• To remove objects with bad photometric data at 



14 



McDonald et al. 



TABLE 5 

Sources showing IR excess in both IOriglia et al.I i|2007I) 

AND this work. 
Td^ Offset^ Temp. Lum. Notes^ 
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LW, A) names from [Clem ent (1997); and 
Lcc — Lcc (197j); 
; van LoonI (12001 ): 



Kg-hand, we remove objects which show an ob- 
served/model flux ratio at Ks-hand which is >0.1 mag 
less than that at both J-band and 77-band (Figure [Ml 
(iv); leaving 127 objects). 

• To remove objects with mid-IR excess caused by miss- 
ing near-IR data (leading to unconstrained fits) , we re- 
move objects with no iiTs-band data unless ([3.6]-[8]) 
> 0.1 mag (Figure UM (v); leaving 120 objects). 

This reduction removes most of the faint giants from 
the sample, but also a few bright giants. These include 
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Fig. 13. — Differences in photometric measurements presented 
in Origlia et al. (2007) and Boyer et al. (2010), showring only those 
stars suggested to harbour circumstellar dust by Origlia et al. 

V19, LWl, and the blended so urces LW7 an d LW8 . Com- 
paring our reduced sample to lOriglia et alj ()2007[ ) shows 
that, out of these 120 sources, 28 match theirs. Of 
these 28, 18 are above 1000 L© an d ten are below 1000 
Lq. The ten objects are listed in lOriglia et al.l (|2007f ) 
as IDs 67, 81, 86, 148 (FBV42), 181, 189 (R7), 222 
(FBV35), 2 9 6, 309 and 373. We note that the study by 
iForte et al.l (|2002D . whose authors lend their names to 
the FBV stars, shows no not ably- different polarisation 



Fundamental parameters and dust production in 47 Tuc 



15 




' Wavelength (urn) 

Fi g. 14. — Examples of stars rejected as having poor photometry. 
See q4.1.2l for details. 



(which would be suggestive of additional, circumstellar, 
line-of-sight dust) in either FBV35 or FBV42, in com- 
parison to the brighter, variable stars. 

As both studies find a significant fraction of potentially 
dust-harbouring stars in the core, we may expect some 
overlap between the two samples, even if they are ran- 
domly distributed. The probability (P) that two samples 
of stars (m and n) are randomly drawn from a total of 
N stars when i stars are found in both samples can be 
found (in the limit of identical object detection efficiency 
and independent datasets) by: 

fc=0 \ \n\{N-n)<. J J 

(2) 

For the entire sample of stars, N — 3399, n — 103 (not 
120, as 17 stars are identified outside the region covered 
by Origlia et al.), m = 93 and fc = 28, giving an answer of 
P <C 10~^°: i.e. zero chance of the overlap being due to 
chance. Concentrating on the cluster's inner 2', however, 
we find P = 0.01% for the bright giantsQ ([8] < 8, i > 
1000 L0) and P = 8.5% for the faint giantfEl (8 < [8] < 
10.874, 35 < i < 1000 Lq). Considering only the inner 
1', P = 0.15% for the bright giantfEl and P = 22.2% for 
the faint giantJ*^. 

We can therefore state that the detection of infrared 
excess around the bright giants has a negligibly-low prob- 
ability of being due to random selection (at a level of 
around 3-4cr). One may also conclude that the detec- 

13 jv = 44, n = 21, m = 22, fc = 16 

14 AT = 254, n = 29, m = 65, k = 10 

15 AT = 24, n = 15, m = 10, k = 9 

16 AT = 131, n = 17, m = 47, fc = 7 




Fig. 15. — Spitzer 8-/^m imaging of the core of 47 Tuc from the 
Rood obs ervations. O range dots mark the stars found to have IR 
excess by lOriglia et al.i 12007 ) ; green circles mark our own dusty 
candidates; large magenta dots denote objects found to have IR 
excess by both parties. Banding is present to the west-north-west 
of bright stars: the brightest, VI and V8, are labelled. 

tion of infrared excess around the faint giants has a non- 
negligible probability of being taken from a random se- 
lection (at 22.2%, or l~2a). However, we remind the 
reader that the two samples are taken from the same 
set of Spitzer data and thus not independent: they will 
therefore have identical noise and image artifacts to con- 
tend with. The real probability that the two samples 
of stars trace objects showing a physical phenomenon 
is (unquantifiably) lower than the 77.8% that the above 
probabilities predict. We are therefore not confident that 
the overlap of 28 stars between our two samples is sig- 
nificant, nor that our two studies are actually tracing a 
particular sub-population of stars. 

4.1.3. Image artifacts 

As noted in M.l.ll the original Spitzer data used both 
in Origlia et al. (2007) and this work contain a number 
of known artifacts: notably banding, the bandwidth ef- 
fect, and blending. Figure [T5l shows the central region of 
the cluster, c o ntain ing almost all the stars suggested by 
lOriglia et al.l (|2007[ ) to contain dust, and also the vast 
majority of our 120 dusty candidates. From this image, 
the lack of overlap between our two candidate sets is 
clear. We now investigate how each of the effects affect 
our data, but stress that we do not discount any individ- 
ual targets in this section. 

Banding is caused by optical scattering and charge 
transfer within the IRAC array, which causes horizon- 
tal and vertical spikes to become manifest around bright 
stars. In Figure [151 these trails can be seen extending 
in the "10-o'clock" and "4-o'clock" directions from the 
brightest stars, with less-prominent lines perpendicular 
to this axis. 

Banding has clearly affected both selections of dusty 
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candidates: several candidates from lOriglia et aLl (|2007[ ) 
lie on banding, caused by VI (IDs 195, 222, 293 and 379) 
amongst others; several of our candidates lie on banding 
caused by V8 (Figure [T5|) . Of the stars both studies 
find to have IR excess, ID-148 and ID-222 are obviously 
affected by banding, though several other stars to the 
west of the cluster may also be affected by this problem 
(IDs 6, 9, 10, 19, 171) as may stars near the center of 
the cluster (including IDs 3, 26, 36, 67 and 181). While 
banding affects stars of all magnitudes, the amplitude 
of the effect scales inversely with received stellar flux. 
Banding can therefore not reproduce the excess in the 
brighter giants, whose excess is also confirmed by either 
their ([3.6]-[8]) colour and/or literature photometry (see 
below, giU). 

Ghost images, caused by the bandwidth effect, affect 
the brightest stars at 8 /im. While this does not appear to 
have affected Origlia et al.'s selection, several of our dusty 
candidates lie on these ghost images, which artificially 
boosts their 8-/im flux, mimicking excess. 

Blending is obviously very significant in the crowded 
cluster core, affecting a set of stars which include the 
aforementioned IDs 3, 26, 36, 67 and 181. Every star is 
blended to some extent, but it is difficult to qua ntify the 
amount of blending present. lOriglia et al.l ()20igi present 
strong evidence from Hubble Space Telescope J-band im- 
agery that very close blends with stars of similar /-band 
magnitude is not an issue in the majority of cases. What 
was not addressed are blends with objects of large (/- 
[3.6]) or (/"[8]) colours, and objects blended with other 
objects at distances of more than 1". 

The majority of very red sources, showing large (/- 
[3.6]) and (/-[8]) colours are galaxies (Paper I). We do 
not, however, consider this a major source of blending in 
this these should affect stars randomly, rather 

than being more prevalent in the inner regions, and are 
likely to only affect one or two s ources in the cluster (cf. 
Paper I; iMatsunaga et al.ll2008[ ). 

On the other hand, blending with objects outside the 
central PSF is still highly significant, as in our exam- 
ple in gTl 2MASS 00234588-7204488. There are two 
main reasons for this. First of all, these images are lim- 
ited by diffraction, rather than seeing, meaning that a 
complex PSF extends a long way from the central source. 
If the flux ratio between the two blended sources is high 
enough, a small error removing flux from the brighter ob- 
ject can leave significant flux in the Airy rings it creates. 
Some of this flux can be attributed to the fainter object, 
causing it to be artiflcially brightened, and leading to an 
obviously non-stellar SED. 

A fundamental difference exists in the reduction pro- 
cess between the Origlia et al. papers and ours. Our 
daomatch/daomaster reduction fits and removes the 
PSF of the brightest sources, then subsequently removes 
fainter objects until we have removed all significant flux 
from the image. Origlia et al.'s ROMAFOT-based reduc- 
tion fits PSFs to the lower-resolution, longer-wavelength 
Spitzer data. They then sum the flux from all sources 
within their K-hand data which lie inside the Airy ra- 
dius of each Spitzer source (2.3" at 8 /im) and use Ku- 
rucz model atmospheres to determine the expected flux 
at longer wavelengths: stars with flux in excess to these 
values are determined to have infrared excess. Each tech- 
nique has its benefits. Our more-conventional technique 



is optimised where the i^-band and Spitzer data have 
similar PSFs (as is the case with the 2MASS data), while 
Origlia et al.'s more-innovative approach should theo- 
retically work better for their higher-resolution /C-band 
data. 

Our concern is whether their method is sophisticated 
enough to take into account all forms of blending: by 
summing flux from sources with centroids within 2.3", a 
blended source at 2.2" from a star is treated the same 
as one lying much closer to the star, while a source at 
2.4" from a star is not considered to affect the star at all. 
In the former case, the blending star will not contribute 
as much to the PSF as is calculated, leading to an over- 
estimation of the expected 8-/xm flux; in the latter case, 
the blending star will contribute some unaccounted flux 
towards the PSF and lead to an under-estimation of the 
expected 8-/xm flux. 

This is of particular concern when a source lies on the 
Airy ring of a bright star, such as Origlia's ID-373, which 
is affected by the nearby VI. One can calculate the con- 
tamination between two stars by convolving the IRAC 
8-/im PSF of the measured star, with the PSF of the 
blending star which has been shifted appropriately in the 
RA-Dec plane and multiplied by the fractional error in 
its flux. One can consider a typical example of two stars, 
where a bright giant is 2.4" from a star 33 x (3.8 mag) 
fainter than it. The bright giant will have a typical error 
in flux of 3%, which corresponds to a flux comparable to 
that received from the fainter star. If one fits a PSF to 
the fainter star, one will find it 0.12-0.14 mag brighter 
than it actually is (depending on the position angle of 
the blended source). It is worth remembering that this 
is caused by a la error: with ^20 giants bright enough 
to cause this effect, one may expect errors up to roughly 
twice this value. When one takes into account blending 
by multiple brighter objects (as happens in the cluster 
core) and any noise inherent in the data, this effect can 
increase further. 

The second problem is that source confusion can lead 
to astrometric errors. A second fundamental difference 
exists between our reduction processes: we flrst flt PSFs 
to objects detected in each band, then combine the data 
from different bands together, while Origlia et al. fit PSFs 
to identical co-ordinates across all bands simultaneously. 
Again, each technique has its benefits: Origlia et al.'s 
method mitigates against displacement of the PSF cen- 
troid due to blended sources, while our technique avoids 
problems with image distortion and alignment within 
surveys, and co-ordinate system alignment among sur- 
veys. 

While the co-ordinates of most of the brightest stars 
agree to within ^0.2" with the positions provided by 
Livia Origlia (private communication), LWIO and V26 
are misplaced by 0.59" and 0.74", respectively. Both 
stars are «2500 Lq, meaning the poor removal of the 
stars' PSFs caused by such an error may lead to signif- 
icant errors in other stars' fluxes. Both studies identify 
several (different) objects around both of these stars as 
having IR excess, but it is not clear that this is actually 
the case. 

Despite the differences between our reduction methods, 
one would expect that if fainter stars had infrared excess 
created by dust, we would mostly identify the same stars 
as having infrared excess. 
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4.1.4. Is {Ks - [8]) better than ([3.6] - [8]) ? 

Of our remaining 120 objects, several exhibit a 'step' 
jump of several tenths of a magnitude between Ks-hand 
and 3.6 /im, but no additional excess at longer wave- 
lengths. Two possibilities for this are large amounts of 
hot circumstellar dust, or bad cross-matching between 
optical/near- IR and Spitzer data due to the very differ- 
ent resolutions between the two. Bad cross-matching is a 
problem that affects all surveys in crowded regions, but 
these objects would be retained under simple colour cuts 
in, say, {Ks-[8]). 

We can separate bad cross-matching from real dust 
by examining where these objects lie in a colour-colour 
diagram (Figure [T6l top- left panel). Strictly speaking, 
this figure is a colour excess - colour excess diagram, and 
thus in the following, we denote a colour excess in filter 
Y with reference filter X as {X-Y)c- This is equal to the 
observed colour minus the model atmosphere colour: 



(X-Y) 



obs 



odd 



(3) 



or, identically, the colour excess of the observation with 
respect to the model in filter X minus the colour excess 
in filter Y: 



observed 



Xr. 



Ddol) 



observed 



Yr, 



jdcl 



)• (4) 



Two populations separate out in this diagram: ones with 
positive colours in both {Ks-[3.6])c and ([3.6]~[8])e and 
ones with positive colours in only (ii's~[3.6])o. 

In Figure [1^] (top-right panel) , we have taken the SED 
of a typical star with {K,-[3.6])c > but ([3.6]-[8])c 
« 0: 2MASS 00234588-7204488. We show circumstel- 
lar dust mo dels made using m etallic iron (optical con- 
stants from lOrdal et al. 119 88*) and amorphous carbon 
(optical constants from IZub ko ct al. 1996): two of the 
best grain species for produc ing flux at short wavelengths 
(e.g. iMcDonald et alJ[20Tol ). Neither can produce the 
observed SED as well as the naked photospheric model, 
multiplied upwards in flux. This is highly suggestive of 
source blending in the Spitzer images, but not in the 
/Cj-band data. Indeed, 'Mc Laughlin et all (j2006f) report 
stars with HST magnitudes of = 14.65 and 15.17 

mag at distances o f 0.335" and 1 - 74" fr om this location, 
respectively, while iSalaris et al.l ()2007D reports the ifg- 
band magnitudes of the same two stars as 12.064 and 
12.018 mag, respectively. On the other hand, 2MASS re- 
ports Kg = 11.485 ± 0.073, whic h compares favorabl y to 
[3.6] = 11.539 ± 0.098 mag from IBover erall (IMoh . It 
would therefore appear that this object is an unresolved 
blend in both 2MASS and Spitzer data. 

A third test we can do to ensure that these stars 
are unresolved blends is to look at the radial distri- 
bution of these objects (Figure \W[ bottom- left panel). 
Sources with significantly positive ([3.6]-[8])o colours ex- 
hibit more central concentration than the main popu- 
lation of bright stars, falling off roughly as the square 
of projected source density, as might be expected from 
blending. Those with near-zero ([3.6]-[8])o colours are 
even more centrally condensed, with the entire group lo- 
cated within 2.5'. 

One can perfo r m a s imilar test with the objects found 
in lOriglia et al.l (|2007[ ). When doing so, we must bear 
in mind the completeness of our samples: comparing the 
distribution of their dusty sources to our total source 



distribution assumes both our datasets achieve the same 
completeness. The fact we work from the same data, 
however, would suggest that their completeness at 8 fim 
is not significantly different to ours. In absolute terms, 
false star tests retrieve 88% of stars with [8] = 10.9 mag 
over the image, though this decreases to 75% of stars be- 
tween 1.5' and 2' and 58% within 1'. These percentages 
increase to 94%, 88% and 87% for stars of [8] = 10.1 
mag. We estimate that we should retrieve ~90% of stars 
of similar magnitude to Origlia et al.'s dusty stars (which 
have an average of [8] = 9.51 mag). Figure [T71 shows that 
the radial distribution of both the bright ( [8] < 8 mag) 
and fain t (8 < [8] < 10.874 mag) giants we detect in the 
original (jOriglia et al.ll2007l ) Spitzer images. 

The distribution of both the bright and faint giants 
follow each other within the Poissonian noi se (-y/rt). The 
distribution of bright dusty candidates in lOrigl ia et al.l 
([2007) matches the global distribution of bright giants we 
identify to within the Poissonian noise. The distribution 
of their faint dusty candidates, however, departs from 
the global distribution of stars with an identical mag- 
nitude range by over 2^yn, showing considerable central 
concentration. The lack of sources at radii > 2' corrobo- 
rates this finding: including stars over 2' from the cluster 
core increases the number of objects by ss41%, suggesting 
there should be ^41 % (27) more f aint d usty sources at 
large radii, whereas lOriglia et al.l ()2007[ ) only find two. 
Similarly, one expects ~15 sources over 2.5' from the 
cluster core, whereas Origlia et al. find none. The dif- 
ferences we find in the radial distribution of their dusty 
objects and our total detected sample (with or without 
correction for unresolved targets) cannot be accounted 
for by differences in detection efficiency. A grossly cen- 
trally concentrated distribution of candidates therefore 
suggests that blending and its ass ociated effect s are a 
significant problem in the sample of lOriglia et al.l (|2007| 1 . 

Finally, we compare the magnitudes fo r all stars 
which are common between the ISalaris et al.l (|2007D and 
2MASS catalogues (Figure [TBI bottom-right panel). The 
majority of both groups of stars (with near-zero and 
with positive ([3.6]-[8])o colours) scatter generally to- 
ward brighter 2MASS magnitudes. The lower-resolution 
of 2MASS therefore suggests that these stars are blends. 
We note that the brightest stars are not in cluded in this 
test, a s they are saturated in the data from ISalaris et al.l 
l|200l . 

To summarise, very warm dust does show an apprecia- 
ble 3.6-/im excess (FigurelTHl). hence a positive (i^s-[3.6]) 
colour. This would suggest that (ifs-[8]) is theoretically 
a better indicator of the presence of warm dust than 
([3.6]-[8]). However, warm dust should still produce a 
substantial ([3.6]-[8]) colour, meaning it should not be 
greatly desensitised to the detection of warm dust when 
compared to (ifs~[3.6]). The above four observations in- 
dicate that the difference in resolution between the Kg 
and [3.6] images has caused artifi cially- high (j^s-[3.6] ) 
colours in both our data and that of lOriglia et al.l ()2010[) , 
while ([3. 6] -[8]) colours remain largely unaffected. We 
therefore conclude that (i^s-[8]) may detect warm cir- 
cumstellar dust more efficiently than ([3. 6] -[8]), but with 
three important caveats: (1) that in dense regions the 
resolution of the images used to determine colour should 
be broadly similar, to avoid unresolved blending affect- 
ing one band; (2) that stars of <3500 K will naturally 
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Fig. 16. — Separation of dusty stars from bad cross-matches and blends. Top-left panel: {Ks — [3.6]) and ([3.6] — [8]) colours for all stars 
in our sample with respect to our model and our colour cut for separating the two populations. 

Top-right panel; SED of source 5.9408750-72.0803889 (black points) which is a poorly-cross-matched source, showing; a model photosphere 
(short-dashed, blue lines) at the fitted luminosity and at 67% brighter than the fitted luminosity, an amorphous carbon dust model at 1500 
K (solid, red line), and a metallic iron dust model at 1500 K (long-dashed, green line). 

Bottom-left panel; radial distribution of (solid, red line) all stars with L > 31.6 with Spitzer photometry, (long-dashed, green line) stars 
with positive colours for both {Ks — [8]) and ([3.6] — [8]), (short-dashed, blue line) stars with positive {Kg — [8]) colour but zero ([3.6] — [8]) 
colour. 

Bottom-right panel; Ks-hand magnitudes of stars covered by both lSalaris et at] II2007I ) and 2MASS. Dusty candidates below the line in the 
upper-left panel are shown as large, black dots, while those above the line are shown as large, blue circles; the line denotes parity between 
surveys. (Note; the brightest stars are not present in the Salaris data due to saturation problems.) 



show a (ii's-lS.G]) colour simply by virtue of being cool; 
and (3) neither colour can mitigate against the effect of 
blending with a red object, such as a background galaxy. 

4.1.5. Examination of individual stars 

With 120 targets remaining, it is now feasible to per- 
form a visual inspection of each SED. We find that many 
faint stars have been selected merely because of their 
(i^<5-[3.6]) colour, which suffers from the problems high- 
lighted in H4. 1.41 and Figure[Tni This point is quite simple 
to address: we apply the colour cut shown in the upper- 
left panel of Figure [111 to exclude stars with ([3.6]- [8] )c 
<0.1 mag which have no 24- /im data. This removes 40 of 
our 120 targets, including Origha et al.'s ID-148, ID-181, 
ID-189, ID-296 and ID-309. 

We now visually identify each of the remaining 80 
targets on the Spitzer 8-/Ltm mosaic and remove those 
objects which lie on bands created by stars, unless the 
photometric excess observed is significantly greater than 
the flux in the band. We also remove targets lying on 



the first Airy ring of bright (>1000 Lq) stars, accord- 
ing to the same criterion, stars artificially brightened by 
bandwidth-effect ghosts, and stars on the edge of the im- 
age. Banding and blending on the Airy disc removes six 
stars due to VI (including ID-222 and ID-373), five stars 
due to V8, two stars due to V4, two stars due to V27, 
one star due to LWl, one star due to LW7 and/or LW8, 
two stars due to LW9, three stars due to LW13 (one of 
which is also on the Airy disc of A19), and two stars 
due to LW18. In the cases of V27 and LW18, the ef- 
fects of banding are compounded by banding from other 
bright stars in the cluster core. Two stars are here re- 
moved as their Spitzer photometry is uncertain, as they 
lie on the edge of the mosaic. Selected objects which are 
actually banding-effect ghost images are also removed 
from around V21, LW3, LW9 and FBV45. This leaves 
50 candidates, with the only faint stars in common with 
lOridia et all (l20Q7l) being ID-67, ID-81 and ID-86. 

Our remaining candidates now split into three main 
groups: 
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Fig. 17. — Radial distribution of bright and faint giants in our 
data, compared to the sample of bright and faint dusty candidates 
in the works of Origlia et al. All data are normalised to 2' and the 
number (A'^) of stars within that radius is shown for each sample. 

1. 29 bright (>1000 Lq) giants that are known 
variables and which (mostly) have known, 
spectroscopically -confirmed circumstellar dust 
()van Loon et al J 120061; ^ebzelter et al.ll2006f )). We 
examine these in ^4.1.61 

2. 9 moderately-bright (250-1000 Lq) giants which 
are listed due to their positive (frs-[3.6])e and/or 
([8]-[24])e colours. We examine these immediately 
after this list. 

3. 14 faint giants (<100 L©) which have Spitzer IRAC 
colours consistent with warm dust. We examine 
these in the remainder of this section. 

All nine moderately-bright giant stars lie in the 
crowded heart of the cluster. Eight of the nine stars 
(excluding F2, below) have not yet been rejected purely 
because their 24-/im photometry still suggested the pres- 
ence of dust: i.e. they had a >10% 24-/im excess. We 
note, however, that the PSFs of these stars are substan- 
tially more blended at 24 fim (FWHM = 6") than at 8 
^m (FWHM = 2"). Seven of the eight stars (excluding 
Fl, below) share all the following criteria: 

• positive {Ks-[3.6])c and/or ([8]-[24])c colours; 

• blended with several other fainter sources within the 
24 ^m PSF; 

• blended with much brighter objects outside the 24 fim 
PSF, including VI, LWIO and LW13 (only six objects, 
but the seventh only shows a 1.6a excess at 24-/im); 

• have (where available) X^-band magnitudes which dif- 
fer between the 2MASS and Salaris surveys, suggesting 
image resolution is affecting magnitude determination; 

• and no other evidence of circumstellar dust. 

Given the problems caused by heavy blending which we 
have already highlighted in this section, wc now remove 
these seven stars, which include ID-67, ID-81 and ID-86. 
This leaves 43 dusty candidates: 29 of these are above 
1000 Lq, and 14 are below. 

None of the 14 faint dus t y can didates we find are in 
common with lOriglia et al.l (|2007l) . We therefore find it 



unlikely that they are dust enshrouded and are likely 
merely artifacts and outliers present in our own photom- 
etry: from our large number of Spitzer sources (4462) we 
may expect outliers up to R:!3.6(t. To be rigorous, and 
because several of these targets lie outside of the field 
of view covered by lOrigha et al.l (2007) , we explore each 
target (F1~F14) individually: 

• (Fl) 00'*23"52.8* -72°04'33", 4305 K, 271 Lq: this 
star was selected purely on the basis of an 11% excess 
at 24 /zm, however the formal error on the flux is 23%, 
meaning this is not a significant detection. 

• (F2) 00''24™05.8" -72°04'44", 4058 K, 260 Lq: sig- 
nificant excess is present at 3.6 and 4.5 /im, but not 
detected at longer wavelengths. It is very close to a 
star four times its brightness, therefore the excess is 
likely a result of blend. 

• (F3) 00''23™47.7" ~72°04'35", 4939 K, 82 Lq: this ob- 
ject shows some flux excess at 5.8 /j,m, strong excess at 
8/im, but is blended with stars of similar brightness at 
these wavelengths. It may also be affected by banding 
from V8, which is expected at this position but not 
visible on the images due to crowding. 

• (F4) 00''25™39.4" -72°08'21", 4445 K, 73 L©: iden- 
tified by an 11% (1.3-1.4(t) excess at both 5.8 and 8 
/im, the low significance of this excess means this star 
is likely a spurious detection. 

• (F5) 00'*23"54.1'' -72°04'17", 5210 K, 71 Lq: this star 
lies on the Airy rings of two slightly brighter stars and 
shows apparent excess at 5.8 and 8 fiia as in F3 above. 
It is also poorly resolved from a nearby, slightly fainter 
companion, and has poor quality J-band photometry 
due to nearby companions. We suspect the apparent 
29% (3(j) excess is due to blending. 

• (F6) 00''23'"59.8" -72°05'01", 5265 K, 59 Lq: this star 
shows increasing excess at 3.6, 4.5 and 8/im. It is not 
detected at 5.8 /im. Poorly resolved from nearby bright 
stars in Spitzer data, and with inconsistent /-band 
photometry, we suggest that this star shows excess due 
to blending. 

• (F7) 00''24™25.9'' -72°19'37", 5204 K, 51 Lq: flagged 
for investigation due to a 25% {3a) excess at 8 /im, 
photometry at other wavelengths shows considerable 
scatter above the reported errors. It may be a field 
star or unusual object, but the excess (if it exists) does 
not appear to be due to dust. 

• (F8) 00''23™39.2" -72°04'04", 5185 K, 50 Lq: this 
star's PSF overlaps with the Airy rings of several 
nearby, unresolved, much brighter objects at 8 /tm. 
It shows an excess at 3.6 and 8 /tm, but a possible 
deficit at 4.5 /tm. We consider the 18% {3a) excess to 
probably be due to blending. 

• (F9) 00''24™07.0" -72°05'46", 4807 K, 47 Lq: we have 
no near-IR photometry for this star, meaning the pho- 
tospheric contribution to the SED is less constrained. 
It was selected on the basis of a 16% excess at 8 /im, 
however this is at a significance of <2a, so is probably 
a spurious detection. 

• (FIO) 00''24™10.0" -72°03'25", 4765 K, 46 Lq: se- 
lected on the basis of a 34% {3a) excess at 8 /im, this 
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star lies on the Airy ring of an unnamed star, which is 
20 X brighter at 8 /xm, but not luminous enough (725 
Lq) to be selected as a potential bright-star blend ear- 
lier in this section. 

• (Fll) 00''24'"11.4" -72°04'04", 4983 K, 44 L©: this 
star was flagged because of its {Ks-[3.6]) and ([5.8]- 
[8]) colours. We note that the same excess persists 
at 3.6, 4.5 and 8 fim, and that the JHKg fluxes are 
deficient compared to our model. The apparent excess 
in this object seems due to the differing resolutions at 
Ks and 3.6 /J,m (see M.f .4p and a spurious 5.8-//m flux. 

• (F12) 00''26'"04.8^ -72°08'16", 5547 K, 42 L©: iden- 
tified due to a 34% {2(j) excess at 8 /um, the excess in 
this object is likely only a spurious detection. 

• (F13) 00''28™35.1" -71°51'09", 40f0 K, 39 Lg: this 
star shows considerable excess at 5.8 and 8 /xm (6(7 
and 4(j, respectively). It is very cool for its luminosity, 
and is 23' from the cluster core. It seems most likely 
that this star is a luminous member of the SMC which 
was marginally brighter than our luminosity cutoff. 

• (F14) 00'*24"3I.7" -72°06'00", 4576 K, 35 Lq: iden- 
tified due to apparent excess at 3.6 and 5.8 /im (the 
only two Spitzer detections), this star has photom- 
etry which departs significantly from the models at 
IJHKs- The SED of this object is of insufficient qual- 
ity and coverage to determine whether any IR excess 
truly exists. 

In summary, of the above 14 stars, we find that the 
excess observed in eight is of low statistical significance 
(< 3.6a), that two have poor-quality SEDs, that three 
suffer from blending and that the final star is likely a 
member of the SMC. We remind the reader that, in any 
case, no ne of these stars wa s found to have mid-IR ex- 
cess bv lOriglia et"all (|20Tol ). To conclude, we find no 
compelling evidence that any star in 47 Tuc produces 
dust until it reaches at least 1000 L©. We find simila r 
mid-IR excesses as those seen bv lOrigha et"all (l200l . 
but we find these in different objects. Such objects oc- 
cur predominantly in heavily-blended environments, and 
are due to either image artifacts or under-reported pho- 
tometric errors due to blending. 

4.1.6. Excess among the bright giants 

We now examine the bright giants which are left as can- 
didates. From these, we remove LW20 and LW17, as they 
were flagged due to their variability, but have no mid-IR 
data with which to determine whether they are dusty. 
Identifying blended objects, we find that LW7 and LW8 
arc blended with each other and that V26 is blended with 
a star just over half its brightness. LW7 and LW8 are both 
identified as having IR excess. It would appear likely that 
at least one of them harbours dust, however the amount 
of excess they show and the quality of their photome- 
try prevents us from saying anything more substantial 
about any such dust. V26 is retained, as the star it 
is blended with is fai nter, and becau se higher-resolution 
8.6-fim photometry (jvan Loon et al.l 2006) confirms the 
excess seen by Spitzer. Five more targets (LW3, FBV45, 
LW18, LW15 and LW6) were removed because their ex- 
cesses are not statistically significant. The remaining 22 
bright giants are listed in Table El 



We again examine the SED and literature data of 
each star in turn. We incl ude in these SEDs litera- 
ture mid-IR pho tometry from Origlia et al.l (I1997D fESO- 
3.6m/ TIMMI); lOridia et al.l (120021 ) (ISO): llta et all 
(|2007[ ) (AKARI) and the AKARI point-source cata- 
logue ()Ishihara et al.|[20Tol ). The ISO and AKARI data 
are taken at coarser resolution than that provided by 
Spitzer IRAC, so caution must be taken when dealing 
with stars near the cluster core. We also include pho - 
tometric data and spectra from Ivan Loon et al.l (|20060 
(ESO-3.6m/T IMMI 2). Comparison with the Spitzer IRS 
data of lLebzelter et al.l (|2006f ) allows us to more-precisely 
determine dust compositions. 

In the absence of spectroscopy, we have seen that iden- 
tification of circumstellar dust production based solely 
on Spitzer photometry can be problematic. Strong IR 
excess which increases in magnitude (but not necessar- 
ily fiux) with wavelength is generally a reliable indicator 
of the presence of circumstellar dust. Spectra are needed 
to conclusively identify the dust composition, though the 
presence or absence of silicates or other oxides can some- 
times be determined in comparatively-isolated stars with 
excess between 8 and 24 /im. Based on the quality of the 
SEDs, literature observations of these stars, and the IR 
excess we observe, we indicate in Table |6] whether a star 
is clearly (Y) or probably (?) dusty. Further analysis 
of the dust minerology and mass-loss rates of these stars 
can be found in the accompanying Paper IV. 

5. DISCUSSION 

5.1. RGB mass loss and metallicity 

Modelling of horizontal branch star masses O implies 
that integrated RGB mass loss per star in 47 Tuc does 
not greatly exceed that in w Cen (see also Paper I) . This 
is despite the clusters being of similar age, and uj Cen be- 
ing (on average) almost a factor of ten more metal-poor. 
Similar HB stellar masses and integrated RGB mass-loss 
rates in the two clusters would imply that the domi- 
nant RGB mass loss process is metallicity- mdependent 
for globular cluster stars. 

This would corroborate findings in other clusters (open 
and globular, including NGC 6791 and uj Cen) that 
the super-solar metallicity open cluster NGC 6791 that 
show RGB and AGB dust production is similar to 
solar-metallicity expectations and that there is no un- 
expe cted absence of RGB stars d ue to super-solar mass 
loss ()van Loon et al.ll2007ll2008af) . This is not to say that 
dust production is necessarily metallicity-independent: 
lower-metallicity stars are warmer at a given luminos- 
ity, therefore their chromospheric mass-loss rates are ex- 
pected to be higher than their higher-metallicity coun- 
terparts (Schroder & : Cuntz 2005 ). Better estimation of 
horizontal branch star masses are required in order to 
improve models before this can be fully substantiated, 
however. 

5.2. RGB dust formation 

We find no reliable evidence that dust is being pro- 
duced by any RGB star below 1000 L©. Th is is co ncur- 
rent with the recent findings of iBover et al.l ()2010D . We 
do not confirm the claims of dust production in lower- 
luminosity RGB stars made in iOriglia et al., (2007., 2010.) . 
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TABLE 6 

Candidate dust-enshrouded stars above 1000 Lq in 47 Tuc. 



RA 


Dec 


Temp. 


Lum. 


TVT 1 

Name 


Dusty? 


Notos^ 


(J2000) 


(J2000) 


(K) 


(Lq) 








00 24 12.64 


-72 06 39.7 


3623 


4824 


VI 


Y 


Silicate dust (see vLMO, LPH) 


00 24 08.57 


-72 03 54.6 


3578 


3583 


V8 


Y 


Silicate dust (see vLMO, LPH); x06 in MvL 


00 24 18.55 


-72 07 58.8 


3738 


3031 


V2 


Y 


Silicate dust (see vLMO, LPH) 


00 25 15.94 


-72 03 54.7 


3153 


2975 


V3 


Y 


Excess in vLMO, LPH; x08 in MvL 


00 24 00.50 


-72 07 26.7 


3521 


2603 


V4 


Y 


Silicate & oxide dust (see LPH), ~RGB-tip 


00 24 07.71 


-72 04 31.4 


3500 


2541 


V26 


Y 


x07 in MvL, =LW13 


00 24 02.46 


-72 05 07.2 


3543 


2324 


LWIO 


Y 


xOl in MvL, also in vLMO 


00 23 50.35 


-72 05 50.3 


3575 


2301 


V21 


Y 


Oxide dusts (see LPH) 


00 23 58.15 


-72 05 49.1 


3374 


2204 


LW9 


Y 




00 24 15.12 


-72 04 36.3 


3510 


2140 


V27 


Y 




00 23 29.98 


-72 22 36.3 


3565 


2122 


Leel424 


Y 




00 24 21.70 


-72 04 13.1 


3526 


2096 


A19 


Y 




nn 0/1 m 


vo no o 
— ( z uo uy.y 


o / io 


zu / y 




V 




00 24 09.39 


-72 04 48.8 


3816 


1640 


MVx3 


Y 


x03 in MvL 


00 24 23.15 


-72 04 22.8 


3738 


1638 


LW19 


7 


x05 in MvL 


00 24 14.44 


-72 05 09.1 


3602 


1575 


V20 


Y 




00 24 08.88 


-72 02 59.5 


3684 


1528 


V22 


7 


JHK photometry questionable 


00 24 25.67 


-72 06 29.9 


3763 


1374 


V6 


7 


No excess in LPH 


00 25 03.60 


-72 09 31.7 


3741 


1363 


V5 


7 


No excess in LPH 


00 25 09.15 


-72 02 39.6 


3692 


1297 


V18 


Y 


Silicate dust (see vLMO, LPH) 


00 24 29.51 


-72 09 07.5 


3775 


1259 


V23 


7 




00 22 58.41 


-72 06 56.1 


3657 


1029 


V13 


Y 


Oxide dusts (see LPH) 



fc WoodI 12005). Other names from 
FBV — IFortc et al.' (WOT); MV —'McDonald & van Loon (2007). 

the software used to determine photometric magnitudes. 
Given the strong artifacts present in the original data 
(Figure [15]), we would advise care in any claim of excess 
among fainter stars in regions affected by these artifacts 
(including the cluster core). Higher- resolution data, such 
as those of Momany et al. (in prep.), are needed before 
the claim of dust production at low luminosities can be 
conclusively proven or refuted, and we await their results 
with anticipation. 

We instead find that dusty mass loss begins in 47 Tuc 
at ~1000 L0 and becomes commonplace by ^2000 L0. 
Between 10(30 and 2000 L©, it is unclear why dust pro- 
duction is not occurring around all stars. In this re- 
gion, nine out of 56 show possible IR excess, while three 
(x03, V18 and V13) show strong mass loss. This could 
represent episodic dust production on either short or 
long timescales, potentially related to pulsation. Alter- 
natively, it could represent a difference in dust produc- 
tion between RGB and AGB stars (due to, say, gravity 
or dredge up), which we cannot separate at these lumi- 
nosities. The luminosity at which dust production first 
occurs in 47 Tuc is therefore very similar to those we find 
in uj Gen (Paper I). 

6. CONCLUSIONS 

In this work, we have used spectral energy distributions 
to establish physical parameters for stars in the globular 
cluster 47 Tuc. We have used these to investigate the 
basic parameters of the cluster as a whole, and the later 
stages of the evolution of its stars. We summarise our 
conclusions as follows: 

• Simple isochrone fits to the cluster's Hertzsprung- 



Lee 



iLia 119771) ; 
&ie the established distance 



and age of the cluster. We find d ~ 4611^200 PC and 
i = 12 ± 1 Gyr. 

• HB models show that mass loss on the RGB is un- 
likely to greatly exceed that in the similarly-aged but 
much more metal-poor cluster w Gen, implying that 
RGB mass loss does not vary with metallicity. We do 
not rule out any correlation of dust production with 
metallicity: our results apply only to integrated mass 
loss. 

• We find that apparent IR excess in stars below L = 
1000 L© is almost certainly due to artifacts in the orig- 
inal Spitzer maps and under-reported photometric er- 
rors ca used by blend i ng, a findi ng contrary to the re - 
sults of lOriglia et"all (|2007[ ) and lOrigha et al.l (|2010D . 

• Some stars above L = 1000 L© show IR excess con- 
sistent with dust production. The fraction of dust- 
producing stars approaches unity above L « 2000 L©. 
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